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Abstract. - We propose a non-BCS mechanism for superconductivity in hole-underdoped 
cuprates based on a gauge approach to the t-J model. The gluing force is an attraction between 
spin vortices centered on the empty sites of two opposite Neel sublattices, leading to pairing of 
charge carriers. In the presence of these pairs, a gauge force coming from the single occupancy 
constraint induces, in turn, the pairing of the spin carriers. The combination of the charge and 
spin pairs gives rise to a finite density of incoherent hole pairs, leading to a finite Nernst signal 
as precursor to superconductivity. The true superconducting transition occurs at an even lower 
temperature, via a 3D XY-type transition. The main features of this non-BCS description of 
superconductivity are consistent with the experimental results in underdoped cuprates, especially 
the contour plot of the Nernst signal. 
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In this Letter we propose a new mechanism of supercon- 
ductivity in hole-underdoped High Tc cuprates, using the 
spin-charge gauge approach to the 2D t-J model, describ- 
ipg the CuO planes [T]. In this approach the disturbance 
of hole doping on the antiferromagnetic (AF) background 
ip systematically considered, giving rise to spin vortices 
dressing the charge excitation (fermionic spinless holon). 
At the same time, due to these vortices the spin excitation 
(■bosonic spin 1/2 spinon) acquires a finite gap, leading to 
4 short range (SR) AF order. The interplay of that SR 
order with the dissipative motion of charge carrier results 
in a metal-insulator crossover, a pronounced phenomenon 
in the underdoped cuprates. A number of peculiar fea- 
tures of cuprates in the normal state can be well explained 
within this scheme [1]. Here this approach is generalized 
to consider the superconducting state. 

The gluing force of the superconducting mechanism is 
an attraction between spin vortices on two opposite Neel 
sublattices, centered around the empty sites (holes), and 
we propose a three-step scenario: At the highest crossover 
temperature, denoted as Tph, a finite density of incoherent 



holon pairs are formed. We propose to identify this tem- 
perature with the experimentally observed (upper) pseu- 
dogap (PG) temperature, where the in-plane resistivity 
deviates from the linear behavior. However, the holon 
pairing alone is not enough for superconductivity to ap- 
pear. Due to the no-double occupation constraint, there 
is a gauge interaction between holon and spinon, through 
which the spin vortex attraction induces in turn the forma- 
tion of spin-singlet (RVB) spinon pairs with a reduction 
of the spinon gap. At the intermediate crossover tem- 
perature, denoted as Tps, a finite density of incoherent 
spinon RVB pairs are formed, which, combined with the 
holon pairs, gives rise to a gas of incoherent preformed hole 
pairs. We propose to identify this temperature with the 
experimental crossover corresponding to the appearance 
of the diamagnetic and Nernst signal. Finally, at an even 
lower temperature, the superconducting transition tem- 
perature Tc, both holon pairs and RVB pairs, hence also 
the hole pairs, become coherent. The proposed supercon- 
ducting mechanism is not of the BCS-type, and it involves 
a gain in kinetic energy (for spinous) coming from the J- 
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term of the spin interactions. The main features of this 
non-BCS description of superconductivity are consistent 
with the experimental results in underdoped cuprates, es- 
pecially the contour plot of the diamagnetic and Nernst 
signal nig. 

The spin-charge gauge approach [T] relies on the follow- 
ing key ideas: 1 ) We decompose the hole operator of the 
t-J model where /i is a spinlcss fcrmionic 

holon, carrying charge, while Za is a spin 1/2 bosonic 
spinon, carrying spin, together with an emergent slave- 
particle gauge field {A^j^), minimally coupled to holon and 
spinon, taking care of the redundant U{1) degrees of free- 
dom coming from the spin-charge decomposition. With 
this choice the no-double occupation constraint is auto- 
matically satisfied via the Pauli principle. 2) lx\ 2D (and 
ID) one can add a "statistical" spin flux (e** ) to ^ and a 
"statistical" charge flux (e**'' ) to h, provided they "com- 
pensate" each other in an appropriate way so that the 
product e~** ft,*e** z is still a fermion. The introduction 
of these fluxes in the lagrangian formalism is materialized 
via the Chern-Simons gauge fields. We find the optimal 
charge and spin statistical fiuxcs in the mean-field approx- 
imation (MFA) [T]. The effect of the optimal spin fiux is 
to attach a spin-vortex to the holon. with opposite chi- 
rality on the two Neel sublatticcs, the rigidity holding up 
a vortex being provided by the AF background. These 
vortices take into account the long-range distortion of the 
AF background caused by the insertion of a dopant hole, 
as first discussed in [S] ; they are naturally associated with 
the semionic nature of the spin fiux (see. [T]) and are rem- 
iniscent of Laughlin's vortices of Ref. [HI Neglecting A- 
fluctuations, the leading terms of the Hamiltonian can be 
written as: 

H = Y.{-t)AM,jh*h,e^^"-^"^ +h.c. 

(iJ) 

+J(1 - h*h, - h*hj)il - |AM,jf ) 
+Jh*h*hjh,\RYB,,\^, (1) 

where AM.y = z|e~'*' e**^ z^- is a kind of Affleck-Marston 
spinon parameter [3] and RVB^ = ^ ^apZiaZjp is an RVB 
spinon singlet order parameter. 3) We use the following 
improved MFA: in the first term of ([T|) we make the MFA 
(AMy ) K, 1, while in the second term we replace the hole 
density by its average, and in the normal state we neglect 
the third term because of being higher order in doping [5). 
Notice that such MFA correctly reproduces the critical 
exponents of the ID t-J model [7], when dimensionally 
reduced (the spin- vortices becomes kink strings in ID). A 
long- wavelength treatment of the second term in ([1]) leads 
to a CP^ spinon nonlinear cr— model with an additional 
term coming from the spin fiux, 

J{V<^'fzh, (2) 

where J = J(l - 2(5) and d^^^'ix) 
e/...9,.^j(-l)|jlA-i(2: - j)h*hj with A the 2D lat- 
tice laplacian. In a quenched treatment of spin 



vortices, self-consistently justified if their density is 
not too small, we derive the MF expectation value 
((V<i>'')^) = mg « 0.5(S| log(5|, which opens a mass gap 
for the spinon, consistent with AF correlation length 
at small 5 extracted from the neutron experiments [S]. 
This gap is also crucial for eliminating an overcounting 
of low-energy degrees of freedom often encountered in 
slave-particle approaches, giving rise to problems in the 
computation of thermodynamic quantities [9]. In fact, 
because of the spinon gap, the low-T thermodynamics 
in this approach is essentially dominated by the gapless 
holons, while the contribution of transverse and scalar 
gauge fluctuations to the free energy almost canceling 
each other |10j . 5^ In the parameter region corresponding 
to the PG "phase" of the cuprates [11] the optimal charge 
fiux is TT per plaquctte and via Hofstadter mechanism it 
converts the spinlcss holons h into Dirac fermions with 
a small Fermi Surface (FS) ep ^ t6. Their dispersion 
is defined in the Magnetic Brillouin Zone (MBZ), which 
we choose as a union of two square regions, denoted 
as R(ight) and L(eft) centered at = (7r/2,7r/2) and 
= (— 7r/2, 7r/2), respectively; these momenta are also 
the centers of the holon FS. Increasing doping or temper- 
ature one reaches the crossover line T* w i/(87r)| log^j 
entering the "strange metal phase" (SM) of the cuprates 
|11| where the optimal charge flux per plaquette is 
instead of tt and we recover a "large" FS for the charge 
excitations with ep ^ t(l — S). 6) Holons and spinous 
are coupled by the dressed gauge field A giving rise to 
metal-insulator crossover and overdamped resonances for 
holes and magnons with strongly T-dependcnt life-time 

m- 

Let us now turn to superconductivity (SC). The SC 
order parameter is assumed to be RVB-like: A'^ = 
^^ap^apCaiCjSj) . lu our MF treatment, neglecting gauge 
fiuctuations it is given by the product of (X]a/3 ^apZaiZpj) 
and {h*h*). Hence both expectation values should be non- 
vanishing to have SC in MFA. 

Holon pairing. In our approach the only quartic term in 
HI) is the RVB-like term which is repulsive. However, spin 
vortices centered on holons have opposite vorticity on the 
two Neel sublatticcs and that produces a long-range at- 
traction, previously neglected in the MFA; this will be our 
key attractive force. Physically it is due to the distortion 
of the AF background caused by the holes. This effect in 
the simplest form was first realized by S. Trugman |12j in 
the early days of High Tc studies, who correctly pointed 
out that putting two holes next to each other on two Neel 
sublatticcs would save energy 2J. We include this effect 
in MFA by introducing a term coming from the average of 
zh in dSl): 

J{z*z) ^(-1)1^1+1^1 A-i(z - j)h*h,h*h,, (3) 

where A is the 2D lattice laplacian. In the static ap- 
proximation for holons ([3]) describes a 2D Coulomb gas 
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with couphng constant g = J{z^z), where (z^z) ^ (A^ + 
mg)^/^ — Too with A « 1 as a UV cutoff, and charges ±1 
depending on the Neel sublattice. For 2D Coulomb gases 
with the above parameters a pairing appears for a tem- 
perature Tph ~ 5/27r, which turns out to be inside the SM 
"phase" [13]. Hence the whole PG "phase" lies below Tph- 
However, we will discuss only the SC arising from the PG 
phase, anticipating that extrapolation to SM phase will 
introduce only quantitative changes. 

To implement the above ideas it is useful to distinguish 
holons on the two Neel sublattices, denoting them as a 
and b, and the two square regions R and L of the MBZ. 
The corresponding holons are called a^, 6^ and a^,b^, re- 
spectively. We also measure the momentum from the left 
and right centers of the two regions Q^, . Since not all 
vortices form pairs, a finite screening effect persists and 
the gas of vortices still have a finite correlation length, 
which we denote as ^ w [Jkp)~'^/'^ [T3]. We keep track of 
the screening effect by replacing (in the long wavelength 
limit) A^"'^ in ([3]) by an effective potential between vortices 
on different Neel sublattices, whose Fourier transform ex- 
tends around in a region of order kp and there it is 
given by Ve}f{q) ~ g/{q^ + ^^^)- The coupling between 
the L and R regions occurs only through Ve//, but since 
VeffiQ^ + Q^) < Veff{0) we neglect this coupling. 

We define: 

A^(p) = j dW,ff{q-p){b%^'i) 

with a = R,L. As in Ref. [151 [16], the d-wave pairing 
symmetry is composed of two p-wave pairing within the 
left and right Dirac cones [T7] in the form 



Ai^flA-} ^''{k){k^-ky), ]ia^ R, , . 

A'*(/c)(-fc^-fcy), ff a = ^' 



We adopt the BCS approximation and the energy spec- 
trum takes the following four-branch form 



S„(fc) = ±\/(Ai±|2ifc|)2-f |A^(fc)|2, a = R,L. (5) 

Neglecting branches not crossing FS and following the pro- 
cedure developed for a spin- wave attraction mechanism in 
Ref. [TMTB] , we obtain for the solution of the gap equation 
on the FS 



A" = A''(fcJ.)«5Cexp{- 



<?^2fc. 



-}, 



(6) 



with a constant c. As common for non- weakly coupled at- 
tractive Fermi systems, the MF temperature at which A'' 
becomes non-vanishing should be identified with the pair- 
ing temperature Tph ■ One reinserts gauge fluctuations and 
recovers gauge invariance in the nodal approximation by 
a standard recipe as discussed in [TB] : One introduces ex- 
plicitly the coordinate-dependent argument of A'*, arg A''. 



The first component of the nodal holon is multiplied by 
g-«2^''g^ , while the second component by e* 2 ^'s A _ r^]^^ 

resulting field is a slave-particle gauge- invariant, hence 
physical, "nodon" . One then reinserts appropriately the 
gauge-invariant vector field = ~ 5^)^ arg A'' by the 
Peierls substitution. It is now easy to derive the low- 
energy effective action in the nodal approximation for a, 
as a variant of QED3. 

spinon RVB pairing. This can be materialized by the 
RVB term in ([T]) combined with the gauge attraction. The 
RVB term is repulsive and in MFA of holons it is irrelevant 
in absence of the holon pairing, but it becomes relevant as 
soon as holon pairing appears. In fact, as shown below, 
the gauge attraction then favors an RVB condensation of 
short-range spinon pairs. Introducing a complex RVB- 
Hubbard-Stratonovich gauge field A^^^^ and treating the 
holon pair in MFA, the RVB-term in ([T]) becomes: 



E 



AS AS* 

(a) in) 

^;^f^^+A^*,e.,z„z,,+/..c. (7) 



In the continuum limit following [18] we present A^^^.> as 
a product of a space- independent, but direction dependent 
modulus factor times a space-dependent phase. 

Neglecting gauge and phase fiuctuations and assuming 
rotational invariance, from ([7]) we derive the modified four- 
branches spinon dispersion: 



uj{k) = ±d{ 



lA-^ 



{J\k\ ± |A« 



(8) 



where to^ — JmQ. The positive branches of the dispersion 
([5]) are similar to those found in a plasma of relativistic 
fermions [19]. It suggests the following interpretation: if 

1 A* I ^0 the spinon system contains a gas of RVB spinon 
pairs, an analogue of Coulomb neutral pairs in the rela- 
tivistic plasma, either in the plasma phase, if (A*)=0 or 
in a condensate, if (A*) ^ 0. For a finite density of spinon 
pairs there are two (positive energy) excitations, with dif- 
ferent energies, but the same spin and momenta. They 
are given, e.g., by creating a spinon up and by destruct- 
ing a spinon down in one of the RVB pairs. Notice that 
the minimum at J\k\ = |A*| in the lower branch is like 
the roton minimum in superfluid helium and has an en- 
ergy lower than m^; it implies a backflow of the gas of 
spinon-pairs dressing the "bare spinon" . Hence RVB con- 
densation would lower the spinon kinetic energy. However, 
to make it occur one needs the gauge contribution. 

For (A*) ^ the global slave particle symmetry is bro- 
ken from J7(l) to Z2, due to the condensation of the charge 

2 e holon and RVB pairs. The Anderson-Higgs mechanism 
then implies a gap '--^ |A*| for the gauge field A^. To ex- 
hibit it explicitly one calculates from ([T]) the gauge effec- 
tive action, obtained by integrating out the spinon. The 
result up to quartic terms for the lagrangian density is 
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with X = ^(argA'' — argA'') a phase field, slave-particle 
gauge- invariant, hence physical, whose gradient describes 
the potential of standard magnetic vortices. From previ- 
ous results one can derive the gap equation for jA'*! in the 
continuum MFA, hence neglecting Xt T = 0, (setting 
A = l), 



1 



_J_ 

dudk- 



(w2 



s\2 



(10) 



In the l.h.s. of Eq. ([TU)) the first term originates from 
the gauge action due to the lowering of the spinon mass 
{iris {m-l — |A*|2)i/2)j while the second term comes 
from the original repulsive Hciscnberg term. The r.h.s. is 
due to spinous; the contribution of the spectrum of gauge 
bosons is neglected as subdominant. Whereas in the slave- 
boson approach the RVB pairs are made of fermions and 
the Heisenbcrg term is attractive, so the pair-formation 
is BCS-like, in our approach the RVB pairs are made of 
bosons, and in our chosen representation the Heisenberg 
term is repulsive, while the pair formation arises from the 
decrease in the free energy of spinous, via the lowering 
of their mass gap, induced by holon-pairing through the 
gauge field. Notice that the leading part of the original 
Heisenberg term is used to provide the AF action for the 
spinous, using the identity (holding for bosonic spinous) 
|^M|2 + \RVB\'^ = 1, see eq. ©. Only the subleading 
term proportional to the holon-pair density is used to ob- 
tain the formation of a finite density of RVB-pairs in (fTO|) . 
so the derived superconductivity appears a little bit in 
the spirit of Laughlin's Gossamer superconductivity [5D] . 
One easily realizes that in (jlOp a non- vanishing solution for 
I A''! is possible only if the second term is sufficiently small, 
i.e., for sufficiently large MF holon pair density. Extension 
to finite T is straightforward including also the contribu- 
tion from the gauge bosons (with an approximate form of 
the effective action due to holons) and roughly estimating 
\{hihj)\ ^ A'^kFJ~^ up to rescaling. The temperature at 
which |A^| becomes non- vanishing is Tps, not yet the true 
condensation temperature, Tc. Notice that from (|T0|) it 
follows necessarily Tp^ < Tph- The d behavior obtained 
from ^ and the finite T version of (|10p for the crossover 
temperatures Tph and Tps, as well as the contour plot for 
different A" are shown in Fig. 1. 

SC transition. The SC transition appears as a XY-type 
transition for magnetic vortices. In fact in the gauged 
XY model (HU]) if the coefficient, jA^p , of the Andcrson- 
Higgs mass term for a is sufficiently small the angular field 
X fluctuates so strongly that it docs not produce a mass 
gap for Op and (e'^) = 0. This is the Coulomb phase of 
the gauged XY (or Stueckclberg) model, where a plasma 
of magnetic vortices-antivortices appears. In the presence 
of a temperature gradient a perpendicular magnetic field 
would induce an unbalance between vortices and antivor- 
tices, giving rise to a Nernst signal. Therefore we con- 
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Fig. 1: The T — 5 phase diagram of the mean field gap equa- 
tion of spinon for different values of MF spinon pairing A° 
(gray lines) which could be compared with different levels of 
the Nernst signal; = is Tps- (The curves at high dopings 
are not quantitatively reliable as they do not take into account 
the crossover to SM). The dashed line is Tph, the "upper PG 
crossover temperature". The dotted line is the crossover tem- 
perature between the PG and SM phases, T* (see also [13p. 
The temperature and are in units of J. 



jecture that this phase of model ((TU)) corresponds to the 
region in the phase diagram characterized by a non-SC 
Nernst signal and a comparison of the experimental phase 
diagram in and the one derived in our model (see 

Fig. 1) supports this idea. 

For a sufficiently large coefficient |A*|2, on the other 
hand, we are in the broken symmetry phase; the fluctu- 
ations of X are exponentially suppressed and (e*'^) ^ 
at T = or there is quasi-condensation at T > and 
the gauge field is gapped. One can prove that, due to 
the fluctuations of the field argA'*, in our approach a 
gapless gauge field is inconsistent with the coherence of 
holon pairs in PG, i.e., coherent holon pairs cannot co- 
exist with incoherent spinon pairs. On the other hand, 
due to the QED-like structure of holons-gauge action, the 
gauge field cannot be gapped by condensation of holon 
pairs alone; only the condensation of RVB spinon pairs at 
the same time can provide a gap to gauge fluctuations. 
Thus as soon as e^^ (quasi-)condenses the same occurs to 
(hihj) so that SC emerges, since the SC order parameter 
is A'^ = A''/{hihj). It follows that T^ < Tps and the 
transition is of XY-type being triggered by the behavior 
of the XY field x- 

A few comments on comparison of the present proposal 
with other models on superconductivity in cuprates are 
in order. It is clear that our proposal differs in an es- 
sential way from the traditional BCS-Eliashberg approach 
|21j . no matter whether the electron-phonon interaction 
or the antifcrromagnetic fluctuations serve as the pairing 
glue. The SC transition occurs here in a similar way as 
in the preformed pairs formalism [22]. From the physical 
point of view, this approach is an implementation of the 
basic idea advocated by P.W. Anderson, attributing su- 
perconductivity to the strong correlation effects in doped 
Mott insulators p5M25j . It shares some similarities with 
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other formahsms exploring the same physical idea, with, 
however, some substantial differences. Both in the stan- 
dard slave-boson [5^ and in the bosonic-RVB phase-string 
[27l[23 approaches the Nernst effect and SC occur due 
to Bose-Einstein condensation (BEC) of bosonic holons. 
Since BEC persists for arbitrary small density in these 
approaches both Nernst effect and SC at T = occur 
as soon as the long-range AFO disappears. The same also 
happens in the standard "preformed pair" approaches [22] , 
due to the persistence of condensation of pairs in the ex- 
treme BEC limit. Instead in our approach the repulsive 
interaction between spinous prevents the appearance of 
the Nernst effect below a critical doping, and the hole 
pairing occurs only when the holon pair density is suffi- 
ciently large to "force" the RVB spinon pairing via gauge 
coupling, while an even higher doping at T = is neces- 
sary to get SC. Similar "critical" dopings also appear in 
the phase- fluctuation approach of (SS], the main physical 
difference with ours being in that approach nodes appear 
in the Nernst phase, whereas in ours a finite FS still per- 
sists and nodes appear only in the SC phase. Furthermore, 
in previous approaches no clear evidence of our additional 
crossover T* appears, as distinct from our Tph (often de- 
noted there as T*). 

Remark: Our approach presents 3 distinct crossover 
lines which different authors have alternatively considered 
as the "pseudogap" crossover: the highest one in T, Tp^ 
where holons start to pair reducing the spectral weight of 
the hole (33) : a second lower one, Tps where incoherent hole 
pairs are formed, mainly affecting the magnetic properties 
since a finite FS still persists; and a third one, T*, inter- 
secting Tps , where one crosses from a large to a small holon 
FS, with a complete suppression of the spectral weight 
for the holes in the antinodal region, with physical effects 
observable experimentally both in transport and thermo- 
dynamics [To]. The first two crossovers have mainly a 
magnetic origin, in the formalism described by the spin 
flux while the third one has a charge origin, due to 
the charge flux and appears only in two-dimensional 
bipartite lattices, see [T]. 

Conclusions: We have proposed a non-pcrturbative 
pairing mechanism in High Tc cuprates. We believe it cap- 
tures the most pronounced characteristics of these com- 
pounds: a strong interplay of antiferromagnetism and su- 
perconductivity. The same Heisenberg interaction derived 
from the strong on-cite Coulomb repulsion is responsible 
for both antiferromagnetism and SC pairing: the lead- 
ing term of that interaction gives rise to antiferromag- 
netism, while its sub-leading term providing the paring 
glue due to vortex-antivortex attraction on the AF back- 
ground. Compared with other proposals on pairing mech- 
anisms, it describes this interplay in a more natural way. 
Although many details of our approach are admittedly 
conjectural, the mechanism of SC proposed here is rather 
complete in its main structure and has the following ap- 
pealing features: I) It is not of simple BCS structure, in 
agreement with some experimental data in underdoped 



cuprates [5nH5^ . 2) SC appears only at a finite doping 
above the critical point where long range AF disappears. 
3) It allows vortices in the normal state, as in the pre- 
formed pair scenario, supporting a Nernst signal. 4) The 
appearance of two positive branches in the spinon dis- 
persion relation for a suitable spinon-antispinon attrac- 
tion induces a similar structure for the magnon dispersion 
around the AF wave vector |33| , reminiscent of the hour- 
glass shape of spectrum found in neutron experiments |34j . 
Furthermore, since the spinon gap has a maximum in J, 
the energy of the magnon resonance is also expected to 
have a maximum. 5) In the SC state the gauge gap de- 
stroys the Reizer singularity which is responsible for the 
anomalous T-dependent life-time of the magnon and elec- 
tron resonances in the normal state. Therefore one expects 
that in the SC state, these resonances become sharper 
at the superconducting transition. The compositcncss of 
excitations within the gauge approach, the holon-spinon 
and spinon-antispinon composites with a gauge glue com- 
ing from the single-occupancy constraint, proved to be 
essential in interpreting the transport [I] and thermody- 
namical [TU] properties of cuprate superconductors, turns 
out to be also crucial for the superconductivity to ac- 
tually occur. Furthermore, the nature of the three-step 
crossovers/transition needed for superconductivity, pro- 
posed in this approach: BCS-like for holon-pair formation, 
kinetic-energy driven for RVB spinon-pair formation and 
XY-like for the true superconducting transition, can be 
seen as specific experimental predictions of this approach. 
A more complete presentation of our approach to super- 
conductivity will appear in |35j . 
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